ANALYTICAL ASPECTS OF THE BROWNIAN 
MOTOR EFFECT IN RANDOMLY FLASHING 

RATCHETS 



DMITRY VOROTNIKOV 

Abstract. The muscle contraction, operation of ATP synthase, 
maintaining the shape of a cell are believed to be secured by mo- 
tor proteins, which can be modelled using the Brownian ratchet 
mechanism. We consider the randomly flashing ratchet model of 
a Brownian motor, where the particles can be in two states, only 
one of which is sensitive the applied spatially periodic potential 
(the mathematical setting is a pair of weakly coupled reaction- 
diffusion and Fokker-Planck equations). We prove that this mech- 
anism indeed generates unidirectional transport by showing that 
the amount of mass in the wells of the potential decreases/increases 
from left to right. The direction of transport is unambiguously de- 
termined by the location of each minimum of the potential with 
respect to the so-called diffusive mean of its adjacent maxima. The 
transport can be generated not only by an asymmetric potential, 
but also by a symmetric potential and asymmetric transition rates, 
and as a consequence of the general result we derive explicit condi- 
tions when the latter happens. When the transitions are localized 
on narrow active sites in the protein conformation space, we find 
a more explicit characterization of the bulk transport direction, 
and infer that some common preconditions of the motor effect are 
redundant. 



1. Introduction 

Brownian motors are nano-scale or molecular machines which can 
produce directed motion when the average force and average temper- 
ature gradient are zero [3l [23]. At first glance, the existence of such 
devices seems to be a paradox, and discretization of the idea really 
did lead to Parrondo's paradox in game theory [21 [7j. Typically, the 
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mechanism involved (ratchet) is based on an interplay between the 
Brownian motion (diffusion), an asymmetric (ratchet-like) potential, 
and nonequilibrium of the system due to chemical or thermal fluctua- 
tions. 

The ratchet principle is ubiquitous and appears everywhere from 
political system to famine cycles, from production strategy to cultural 
studies. Motor proteins, which provide muscle contraction (myosin, 
kinesin, dynein), ATP synthase, as well as membrane-bound motor 
proteins maintaining the shape of a cell, can be modelled using the 
Brownian ratchet mechanism [Tj [TTJJ [HJ ETJ [22] • A ribosome can also 
be considered Brownian ratchet device EH1 . 

Motor proteins can attach to and detach from a substrate of vectorial 
symmetry [fl [22] under the action of a chemical energy source. This 
leads to the following boundary value problem [U [TOj [161 fTT] : 



(1) 



Pt ~ °~Pxx - K>(lf) x p) x = vP 
Pt - SPxx = ~VP + VP, 

ap x + Kif> x p = 0, 
P x = 0, 
p>0, P > 0, 



rjp, xe (0,1), t>0, 

x e (o,i), t > o, 

x = 0, 1, t > 0, 
x = 0, 1, t > 0, 



/ p(t, x) + P(t, x) tte = 1, t> 0. 



Here x) and P(t, x) are the unknown densities of the particles in 
"attached" and "detached" states, resp., at a time t and a spatial point 
x; a and q are the diffusion coefficients of "attached" and "detached" 
particles, resp.; ip(x) is the potential; k is a coefficient inversely pro- 
portional to temperature; v(x) > and rj(x) > are the rates of 
transition from one state to another, i.e. v indicates the probability of 
seizing a motor protein by a "detached" particle (located at a spatial 
point x), and rj expresses the probability of losing its motor protein 
for an "attached" particle. A typical "ratchet-like" potential ip with k 
teeth, k > 1, is 1/fc-periodic in x and has a unique local (and, hence, 
global) minimum within each period. 

Let us compare this model with a somehow simpler one called the 
flashing ratchet [3]. Here there is only one state, p is the unknown 
density of particles, a is the diffusion coefficient, and if) is the potential, 
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which is switched on and off cyclically: 



(2) 



Pt = crpxx + h(t)(if) x p) x , x G (0, 1), 
ap x + h{t)ip x p = 0, x = 0, 1, 
P>0, 



f p(x, t)dx = 1, 
o 

h(t) = 1, nT < t < nT + T tr , n = 
h(t) = 0, nT + T tr <t<nT + T, 



0,1,.. 

n = 0, 



Here each particle is potential-sensitive and potential-insensitive for a 
priori known moments of time, whereas in the first model this is deter- 
mined by random attachments and detachments of the motor protein. 
Therefore model ([1]) is sometimes referred to as the randomly flashing 
ratchet. 

Another model related to ([I]) is the collaborative ratchet [5J QUI III 
2Ij [22] where the particles in two states are sensitive to two different 
potentials ip and \1/ which help each other to achieve the motor effect: 



(3) 



'Pt ~ o-p xx - n{ip x p) x 
P t - qP xx - K (V x P) a 
ap x + nip x p = 0, 
sP x + k^ x P = 0, 
P > 0, P > 0, 



vP — rjp, 

= —vP + Tjp, 



x e (0,1) 

x g (o,i; 

x = 0, 1, 
x = 0,1, 



,t>0, 
,t>0, 
t > 0, 
t > 0, 



fp{t,x) 
o 



P(t,x)dx = 1, t > 0. 



This model is relevant not only in connection with biology, but also in 
transport of cold rubidium atoms jH]. 

A diffusion-free problem somehow similar to (jSJ) is a model of trans- 
port of neurofilaments in axons; its behaviour near equilibrium was 
investigated in [6]. 

The mathematical studies of Brownian ratchet models start with the 
question whether they indeed generate unidirectional transport (which 
is observed in experiments and simulations), and what assumptions 
are needed for that. For example, model ([1]) is too general to produce 
transport with any choice of parameters: setting 

rj/u = exp (nip/a), 

we observe that the total amount of mass in the "wells" of the potential, 
i.e. in the segments [|, i — 0, . . . , k — 1, eventually with the course 
of time tends to 1/k. This shows that the ratchets should be "tuned" 
to work well. 
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Comparing eventual distribution of mass between the "wells" of the 
potential -0, it is possible to mathematically endorse the motor effect 
[12] . A left-to-right chain of inequalities in this distribution would mean 
unidirectional transport. 

The occurrence of transport in model ([3]) for small a = q and certain 
interplay between the potentials was shown in [J]. Some of the results 
of that paper — in particular, existence of unique solutions to ([3]) and 
to the stationary problem 

-crpxx - k(iI>xP)x = vP -T]p, x e (0, 1), 

-qP xx - k(V x P) x = -vP + VP, are (0,1), 

ap x + Ki[) x p = 0, x = 0, 1, 

qP x + k^ x P = 0, x = 0, 1, 

P > 0, P > 0, 
i 

/ p(x) + P{x) dx = 1, 



(4) 



and eventual convergence of solutions to (j3J) to the solutions of (jlj) - 
are valid for ^ = const as well. The transport result was generalized 
to the multi-state systems with several interacting non-flat potentials 
in PES]. 

The analytical proof of the motor effect for the flashing ratchet (J2j) 
was given in [25], based on a framework developed in [12]. The uni- 
directional transport occurs when the potential is asymmetric, and its 
direction is determined by the location of the minima of the potential 
with respect to the centres of the corresponding wells. A homogeniza- 
tion approach to the flashing ratchet (j2J) was proposed in [2U]. This 
approach was applied to the randomly flashing ratchet flTJ in [TS] (see 
also [15J). There were presented examples of transport in the case 
of a sawtooth potential (which is asymmetric) and constant transition 
rates (which are obviously symmetric), and with a symmetric potential 
and asymmetric transition rates. We are not aware of any works with 
rigorous mathematical evidences of transport for the non-homogenized 
problem flTJ). 

In this paper, we show that the motor effect in model ([T]) is due 
to a subtle interplay between the asymmetries of the potential ip and 
the transition rate v. More precisely, an asymmetry of the transition 
function v yields a deviation of the so-called diffusive means of the edge 
points of the wells of the potential — these edge points are the maxima 
of ip — from the centres of the wells. The direction of transport is 
determined by the location of the minima of the potential with respect 
to that biased centres (the diffusive means of the adjacent maxima). 
The ratchet is tuned when the influence of the potential on the particles 
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which carry the motor protein dominates their diffusion (which is a 
natural assumption since these complex particles are larger and thus 
diffuse slower). 

Among the consequences of the general result, the following obser- 
vation is of special importance in connection with motor proteins. The 
function v(x) is conformation-dependent, and is frequently assumed 
[TO] to dramatically increase on narrow regions corresponding to the 
areas in the conformation space where the reattachment activity of 
motor proteins is overwhelming, and to be negligibly small everywhere 
else. In this particular framework, we infer that the bulk transport di- 
rection is dictated by the location of the minima of the potential with 
respect to that active sites, and observe that several preconditions of 
the motor effect found in biophysical literature are unnecessary. 

The paper is organized as follows. In Section 2, we set the problem 
more rigorously, define the notion of the diffusive mean, and formulate 
the main result (Theorem 12. ip . Section 3 is a discussion of conse- 
quences of Theorem 12.11 In particular, we investigate the case men- 
tioned in the previous paragraph. In this section, we also study some 
properties of the diffusive mean and evaluate it. In Section 4, we intro- 
duce a semidiscretized device, the squeezing ratchet, and prove that it 
generates unidirectional transport. The proof of the main theorem is 
provided in Section 5, where we show that the squeezing ratchet and the 
original randomly flashing ratchet have similar behaviour. Section 6 is 
a biologically-oriented discussion of the results. We opted for moving 
all the proofs except the one of Theorem 12.11 to several appendices. 



2. Preliminaries 



We consider the stationary boundary value problem for the randomly 
flashing ratchet equation with Neumann boundary conditions 



(5) 



-qP xx = -VP + T]p, 

ap x + ni) x p = 0, 

Px = 0, 

P > o, P > 0, 
1 

J p(x) + P{x) dx = 1, 



VP, 



xe (0,1), 
xe(0,i), 

x = 0, 1, 
x = 0, 1, 



which describes the eventual distribution of particles subjected to the 
action of the ratchet. 

The potential ip(x) and the transition rates v{x) > and t](x) > 
are assumed to be smooth scalar functions on [0, 1] of period 1/k, with 
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k > 1 being a fixed integer. The potential ip should have maxima at 
points Xi and minima at points a.j, and be monotonic (without zero 
slopes) between these points, where 

i — 1 

(6) Xi = ——, i = l,...,k + l, 

k 

(7) ^ — a + Xi, i — 1, . . . , k. 
The positive parameter a should be less than 1/k. 

B 

We use the bra-ket notation (/x, /) = / f dfi, where \i G C*[A,B], 

A 

and / is a continuous function on [A, B]. 

The symbol 5 X denotes the Dirac delta centered at x G R. The sym- 
bol C will stand for a generic positive constant that can take different 
values in different lines. 

Definition 2.1. Unless otherwise specified, we say that u e L\(A,B) 
is a solution to the problem 

f {k\u) xx (x) + (A^zOe) + h(x)u(x) = Q(u)(x), x G (A, B), 
\ (hu) x (A) = (kiu) x (B) = 0, 

where functions ki,k,2, ks G C[A, B], (vl) = ^(-B) = 0, and a linear 
operator 6 : Li(A,B) — > C*[A, B] are prescribed, if 

B 

(8) y — ki{x)u{x)ip x {x) + A;3(x)m(x)<^(x) (ix 

A 

= (e(u),<p) 

for any G C 2 [A, 5], = <p x (B) = 0. 

Note that OH]) already includes the Neumann boundary condition. 

Definition 2.2. Let A and 5 be real numbers, and 4>{x) > be a 
continuous scalar function on [A, B\. The number s G (A, B) is called 
the 4>- diffusive mean of A and B provided the system 

( <pU -U xx = 6 B , x G (A, B), 

(9) <^ C/(A) = [/(S), 

[ U X (A) = U X (B) = 

has a solution C/. 

Proposition 2.1. TTie 4>-diffusive mean always exists and is unique. 
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For any integrable scalar function <3> on (0, 1), we denote 



(10) 



3>(x) dx, i = 1 k. 



The main result of the paper is 

Theorem 2.1. Let S be the u/q- diffusive mean ofx\ = andx2 
If a < S, then, for sufficiently small a and sufficiently large n, 

(11) Pi > P2 > ■ • ■ > Pk, 

(12) P-y > P 2 > ■ ■ ■ > P k . 

Remark 2.1. We will even prove that 

(13) P{x) > P{x + l/k), < x < 1- 1/Jfe, 
which is stronger than (|T2l) . 

We can renormalize to get 

-ap xx - n(i/; x p) x = vP - rjp, x G (0, 1), 



(14) 



-qPxx = -VP + Tip, 

ap x + nip x p = 0, 
P x = 0, 

P > o, P > o, 

1 



x G (0,1) 
x = 0, 1, 
x = 0, 1, 



J T]{x)p{x) + u{x)P{x) dx = 2. 
V, o 

Integration of f|14|) implies 

i 

(15) y r/(x)p(x) — u(x)P(x) dx = 0. 

o 

Since and 1 are maxima of i/j, we have ^(0) = ^(l) = 0. Thus, 
is equivalent to 

-ap xx - K(i) x p) x = vP -rjp, x G (0, 1), 



(16) 



-qP xx = -vP + rjp, 
Px = 0, 
Px = 0, 

P > o, P > 0, 
1 

J rj(x)p(x) dx = 1, 



1 

J u(x)P(x) dx = 1. 



xe (0,1), 

x = 0, 1, 
x = 0, 1, 
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Since problem (jSJ) is linear, it is enough to prove Theorem 12. II for the 
renormalized problem (fTTJj) . 

3. Properties of diffusive means and applications of 

Theorem 2.1 

Theorem 12.11 implies that, if the diffusion of potential-sensitive par- 
ticles is slow, and the influence of the potential is strong (or the tem- 
perature is low), and if a < S, then, given any initial distribution of 
density, the mass of particles of each kind in the wells will eventually 
decrease from left to right, i.e. the motor effect is present. In order to 
present more specific consequences of Theorem 12.11 we need to better 
understand the notion of the diffusive mean. 

Example 3.1. Let 

4>(x) = (f)(A + B - x). 

Let s be the 0-diffusive mean of A and B, and let s be the 0-diffusive 
mean of A and B. Then 

(17) s = A + B-s. 

In particular, if is a constant function or merely 

(18) <j>(x) = $(x), A<x<^^-, 

then, due to (|17[) . the diffusive mean s coincides with the arithmetic 
mean ^y^. In other words, s is located in the centre of the segment 
[A,B). 

In general position, the diffusive mean is distorted from the centre 
of the segment. We now evaluate the direction of the bias. 

Proposition 3.1. Let 

(19) <j>{x) < 4>(x), A < x < ^-t^, 

and the inequality is strict at least at one point. Then the <p- diffusive 
mean of A and B is strictly larger than the arithmetic mean . 

A symmetry argument shows that if 

(20) 0(x)>0(x), A<x<^^-, 

and the inequality is strict at least at one point, then the diffusive mean 
is strictly less than ^y^. 

Theorem 12.11 indicates that the direction of transport generated by 
the ratchet is merely determined by the location of a with respect to 
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S (in other words, by the direction of the vector Sa). Indeed, we have 
just pointed that out for a < S, and if a > S, then p% < p 2 < ■ ■ ■ < pk, 
Pi < P2 < ■ ■ ■ < Pk — to see this it suffices to make the change of 
variables x — Y 1 — x, to apply the reasoning of Example 13 . 1 1 with <fi = ^, 

A = and B = 1/k, and to take into account that </>(x) = (f>(% — x) = 
0(1 — x) due to periodicity. 

Example 3.2. Let the potential if) be symmetric, i.e. ip(x) = if>(l — x), 
which can only happen when a = ^r. The one-state flashing ratchet 
cannot generate transport in this case. However, the randomly flashing 
ratchet can. Assume that 

(21) v{x) < i/(-x + l/k), <x < — , 

with strict inequality at least at one point. Then a = < S by 
Proposition 13.11 and the transport (to the left) occurs. If inequality 
( ]2~T]) has opposite sign, the transport is to the right. 

Example 3.3. Let v be symmetric, i.e. v{x) = v{\ — x). Then peri- 
odicity of v and Example 13.11 imply that S — ^r. The motor effect is 
provided by the condition a ^ ^c, i.e. the potential should be asym- 
metric, and the direction of transport is determined by the location of 
a with respect to as for the one-state flashing ratchet [2"B] . 

The rate of attachment v of free motor proteins to the cargo is deter- 
mined by the coupling between the conformation and the underlying 
chemical reactions. A representative situation is when the function 
v{x) is peaked on narrow regions corresponding to the active sites in 
the conformation space, and almost vanishes outside of that regions 
[TU] p. 55]. Due to periodicity, in each well there is one 

active region. The following proposition makes very transparent the 
application of Theorem 12.11 to this framework. 

Proposition 3.2. Let <p n : [A, B] — > (0, +00) be a sequence of contin- 
uous functions. Let s n be the corresponding <f> n - diffusive means of A 
and B. Assume that there exist numbers s* G (A, B) and K* > such 
that 4> n — > K*5 Sii weakly-* in C*[A, B] as n — > +00. Then s n — > s*. 

Consequently, if the active sites are located around the points Sj = 
s*+Xi, i = 1, . . . , k, then the ^/^-diffusive mean of x\ and x 2 is approx- 
imately equal to s*. Thus, the direction of transport is determined by 
the sense of the vector ~s~^h. Hence, the factor which plays the crucial 
role in the transport is not the asymmetry of the potential but the lo- 
cation of the minima of the potential outside of the active regions. To 
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avoid reader's confusion, we again point out that here we are referring 
only to the active sites where the attachment of the motors occurs. 

However, the detachment of the motors may also be localized within 
small active zones. Normally [TJ [10], it is assumed that the areas of 
detachment activity are located around the minima of the potential. 
Moreover, it is generally supposed that the asymmetry of the potential 
should correlate with the transport, i.e. the position of the minimum of 
the potential should be biased from the center of each well in accordance 
with the intended bulk movement direction, as in Example 13.31 We 
have just seen that all these hypotheses are redundant (provided the 
ratchet is well-tuned in the sense of Theorem 12. ip . 



4. Squeezing ratchet 

We divide the proof of Theorem 12.11 into two main steps. Firstly, 
in this section, we consider a somehow limiting case of the randomly 
flashing ratchet. It appears as a pair of purely mathematical state- 
ments, which, however, have a nice physical interpretation. Namely, 
we will propose an idealized ratchet-like engine, and the just mentioned 
mathematical facts will entail that the direction of transport provided 
by this machine is determined by the location of a with respect to S. 
Secondly, in the next section, we will perform some passages to the 
limit which will imply the claimed behaviour of the randomly flashing 
ratchet. 

Let G(x, y) be Green's function of the Sturm-Liouville operator 

C-- — 

dx 2 

on (0, 1) with homogeneous Neumann boundary condition. Let g»(ar) = 
G(x, aj), i.e. 

/ 99 x / m - sq ixx = S atJ x e (o, l), 

{ ] I «J0) = q ix (l) = 0, 

and let 

k 

(23) ? = X> 

i=l 

As we will observe in the proof of Proposition 12.11 (Appendix A), 

(24) q t (x) >0, xe [0,1]. 



The next theorem is our core tool for catching the motor effect. 
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Theorem 4.1. There exists a unique function Q G C[0, 1] solving the 
following problem: 



(25) 



k . . 

vQ-<&xx = l£(yQ)i8ai, ace (0,1) 

Q»(o) = Q«(i)~= o, 
i 

/ u(x)Q(x) dx = 1. 

k 



Moreover, if 

(26) g(» > g(x + 1/ife) + 7, < x < 1 - 1/k, 
with some 7 > 0, i/ien 

(27) Q(x) > Q(x + 1/k) + M7, < x < 1 - 1/Jfe, 
where M = min i= i v .. jfe (vqi) k . 

Remark 4.1. Theorem 14. II may be considered as a continuous version 
of a purely algebraic fact, [251 Lemma 3.2]. 

Our "left drift" condition a < S arises from the following lemma. 

Lemma 4.1. If a < S , then there is 7 > such that fT26|) holds true. 



The results of this section can be interpreted as follows. Consider 
a semidiscretized device which we refer to as the "squeezing ratchet" , 
and which acts as follows. The particles can be in two states, ground 
and excited, and u(x) and rj(x) are probabilities of transition from the 
first to the second state and back, resp. The particles in the ground 
state diffuse with diffusion coefficient q. If a particle positioned at the 
segment (xi,Xi+i) suddenly changes its state from the ground to the 
excited one, then it instantly jumps to the point a« (located to the left 
from the v/q- diffusive mean of the points Xi and £4+1). Then, given any 
initial allocation of particles, the renormalized eventual distribution Q 
of ground particles satisfies (127)) . i.e. their mass is transported to the 
left. Moreover, the excited particles are eventually concentrated at the 

points cii, and one can observe that the asymptotic amounts Xi ~ ( u Q)i 
of excited particles at the points a« decrease from left to right. A 
reflection argument shows that if a« are located to the right from the 
iz/^-diffusive means of the corresponding endpoints X{ and Xj+i, then 
both ground and excited mass is transported to the right. 
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5. ASYMPTOTICS OF THE TIME-DISCRETIZED FOKKER-PLANCK 
EQUATION AND BEHAVIOUR OF THE RANDOMLY FLASHING 

RATCHET 

This section contains a series of limiting procedures which, together 
with the results of the previous section, will lead to derivation of The- 
orem [2]TJ The reader not interested in the details of the proof of the 
main theorem may skip this section. 

Denote by d the Wasserstein metric of order two on the space of prob- 
ability measures on [0, 1], see e.g. [SJ. The convergence in Wasserstein 
metric is equivalent to the weak-* convergence of probability measures: 

(28) d{fi n , fi) & (fi n - /I, f) ->■ o, fe C[0, 1]. 

Set b(x) = ip x (x)/r)(x). Note that b(x) is zero at the extrema a« 
and Xi of the potential ip, is negative for xi < x < a,i, and is positive 
for cii < x < x i+ i, i = 1, . . . , k. We also recall that k is a physical 
parameter inversely proportional to temperature. 

Each of the following results will be used for the proof of the subse- 
quent ones. We begin with a density-based lemma. 

Lemma 5.1. If u G (7* [0,1] satisfies 

(29) (u, if + Khp x ) = 

for any <p G C 2 [0, 1], y? x (0) = ^(1) = 0, then u — 0. 
We now study the asymptotic behaviour of the system 



(30) 
where 



r - K{br) x = R, x G (0, 1) 
r > 0, 



R G Li(0, 1), R > 0, J R{x) dx = l 



i 

is prescribed. Clearly, J r(x) dx = 1, so r can be considered as a 

o 

probability measure. 
Lemma 5.2. We have 

k 

(31) lim sup d(r,y^ Ri8 a .) = 0. 

K— S- + 00 j ' * 

fle£i(0,l), H>0, fR(x) dx=l 1=1 
o 
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We now need to consider the time-discretized Fokker-Planck prob- 
lem: 

w a>K - a{w^ R /rf) xx - K(bw a:K ) x = R a>K , x E (0, 1), 
(32) { (w a Jr)) x (0) = (w a J V ) x (l) = 0, 



W„ K > 0, 



where 



R a>K E Li(0, 1), R a , K > 0, J R a , K (x) dx = 1, 



depends on a and k. We want to study the asymptotics of this problem 
for large k and small a. 

Lemma 5.3. For each k there exists e K > so that 

k 

(33) lim d(wv K ,) (R BrK ).8 ai ) = ®. 

K— ¥+oo, cr<e K ' ' 

1=1 

We are in a position to obtain the main result of this section, namely, 
to prove that the randomly flashing ratchet behaves similarly to the 
squeezing ratchet. 

Lemma 5.4. For each k there exists t K > so that the corresponding 
solutions of f fT6|) have the following properties: 

(34) lim sup \P(x) - Q{x)\ = 0, 

n^+oo, a<e K o<a;<l 

(35) lim ^ (Z7g)i 



k-5>+oo, a<e K V{ a ) 

The proof of Theorem 12.11 is now straightforward. Inequality (J27j) 
yields 

(36) ^>% + C,i = l,. 

77(a) 71(a) 

Therefore, (TTTT) and (fl3j) are direct consequences of Theorem 14.11 and 
of Lemmas 14.11 and 15.41 

6. Discussion 

We investigate the randomly flashing ratchet model for motor pro- 
teins. In this model, every Brownian particle can be in two states: with 
or without a motor protein. The particles carrying a motor protein are 
sensitive to a chemically-induced periodic potential. The motor-free 
particles diffuse normally. The switch between the two states happens 
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when a particle loses or seizes a motor protein, and the probabilities of 
these events are prescribed. 

The real nature of transport generated by motor proteins is much 
more complex, but the model under consideration is flexible enough to 
encompass the four main factors involved: chemical potential, confor- 
mation, attachment/detachment of motor proteins, and diffusion. 

We have rigorously shown that unidirectional transport of mass oc- 
curs when the diffusion of the potential-sensitive particles is strongly 
dominated by the influence of the potential. The direction of trans- 
port is unambiguously determined by a certain interrelation between 
the asymmetries of the potential and of the function describing the 
probability of attachment of motor protein molecules to the cargo (in 
other words, of the rate of transition from the potential-insensitive to 
the potential-sensitive state): in particular, at least one of them should 
be asymmetric to secure the transport effect. 

In the literature it is often assumed [10] that the attachment of mo- 
tor proteins to transportable objects predominantly happens within a 
narrow active region determined by the conformation. This region is 
unique in every period of the potential. As an application of our general 
result, we show that in this specific framework the ratchet mechanism 
moves the transported matter to the left /right provided the minimum 
of the potential within each period is located to the left/right of the 
attachment zone. We do not require additional hypotheses imposed 
in earlier biophysical papers such that correlation of the asymmetry of 
the potential with the intended direction of transport or localization of 
the detachment areas around the minima of the potential. 

Our theory is also able to detect the transport direction in the case 
when there are two or more attachment sites within each period or 
their spatial distribution becomes more homogeneous. 

Another popular model for molecular motors is the classical or de- 
terministic flashing ratchet. In the bio motor context, it may be inter- 
preted assuming that during some period of time all the particles are 
bound to motor proteins, and are therefore potential-sensitive. Then, 
suddenly, all the tiny "parcels" simultaneously lose their carriers and 
start to diffuse freely. After a certain time interval, the motors syn- 
chronously reattach to their cargoes, and the whole process repeats on 
and on. Thus, conformation is completely disregarded. The direction 
of transport is dictated by the location of the minima of the asymmetric 
periodic chemical potential with respect to the centres of the potential 
wells [25] . The deterministic flashing ratchet has promising applica- 
tions in many areas of science. However, in connection with motor 
proteins, the random ratchet seems to be more pliable and realistic. 
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Our analysis shows that predictions of the random and deterministic 
flashing ratchets may significantly differ. For instance, if the region 
of attachment activity is unique (up to spatial periodicity), and the 
minimum of the chemical potential is located between that active zone 
and the centre of the potential well, then the bulk transport directions 
anticipated by the two models are exactly the opposite of one another. 
On the other hand, when the distribution of attachment sites is sym- 
metric (no matter how many active areas are present in each well) , the 
forecasts of the models essentially coincide. These observations give a 
new insight for experimental verification and validation of the flashing 
ratchets. 

Appendix A. Proof of Propositions 12.11 and 13.11 

Let G(x, y) be Green's function of the Sturm-Liouville operator 

C--— + 

dx 2 

on (A, B) with homogeneous Neumann boundary condition. Then 
U(x) = G(x, s) is a solution to dHJ) if an only if G(A, s) = G(B, s). 

By the distributional maximum principle [13, Theorem B], G(x, y) > 
0. Observe that 

X 

(37) G x (x,A) = J </>(z)G(A,z)dz, 

B 



(38) G x (x,B) = J <j>{z)G{B,z)dz. 

A 

Hence, the function G(x, A) is decreasing in x, and G(x, B) is increas- 
ing. Thus, the function 

g(x) = G(A,x) - G(B,x) 

is also (strictly) decreasing. At the ends of the segment, we have g(A) = 
G(A, A) - G(B, A) > and g(B) = G(A, B) - G(B, B) < 0. Since g is 
a continuous function, there is unique s G (A, B) such that G(A, s) = 
G(B, s). The first proposition is proven. 
Now, let 

gi(x) = G( , x) - G( , A + B-x). 



We claim that 

(39) 9l (A) > 0. 
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Since G is Green's function, 

(40) ( 9l UA) = 0, 
and 

(41) (gi) xx (x) 

= <f>( x )G(^-?-, x )-<f>{A + B- X )G{^^-, A + B-x) 

A + B 

= <j){x)g x {x) + [4>(x) - <j>{ x )]G{—^—, A + B-x) 

, , . . . , A + B 

< <P{x)gi{x), A < x < — - — , 

and the inequality is strict at least at one point. In particular, gi 
cannot be identically zero. 

Assume that gi{A) < 0. By the maximum principle, gi cannot have 
non-positive minima within (A, ^y^). But 

(42) 9i(^-) = 0, 

so A must be a minimum point. Let #2 (a?) = gi(x) — gi(A). Then g2 is 
non-negative, and 

A + B 

(43) (92)xx(x) < <j){x)g 2 {x), A < x < 
Thus, 

(44) (g2) x {x) < [ <j>(t)g 2 (t) dt, A<x< 



Since g 2 {A) = 0, by the mean value theorem, 
(45) g 2 (x) = g 2 (x) - g 2 (A) = (x - A){g 2 ) x {c) 

C X 

<{x — A) J (j>{t)g 2 (t) dt< l -{B-A) J <f>{t)g 2 {t) dt 

A A 

for some c, A < c < x < ^y^. The Gronwall lemma implies g 2 = 0, 
so, by (142|) . gi = 0, and we get a contradiction. 
Hence, 

9(^-) = 9i(A) > 0. 

We have observed above that g(B) < 0, so there is s G (^^-,B) such 
that g(s) = 0, and this number s is the 0-diffusive mean. 
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Appendix B. Proof of Proposition 13.21 



Without loss of generality, [A, B] = [0, 1]. Denoting the solutions of 
)]) corresponding to (fi n and s n by U n , we infer that 
i 



(46) 



( j )ri {x)U 2 n {x) + {U n )ldx = U ri 



The function U n is positive and thus convex on (0, s n ) and (s n ,l). 
Hence, due to the boundary conditions, the derivative (U n ) x (x) is pos- 
itive/negative when x < s n / x > s n , and tends toward its supre- 
mum/infimum as x approaches s n from the left/right, resp. Moreover, 
the maximum of U n is achieved at s n , and the minima are reached at 
and 1. 

It is easy to see that 

i 

(47) min U 2 n (x) < / n (x)t/ n 2 (x) dx 

xe[o,i] A, J 



for large n, and 



(48) 



max Ut (x) < 2 
xe[o,i] ~ 



min U„ (x) + 

xe[o,i] 



(U n )l dx 



Hence, (JIB]) implies 

(49) 
whence 
(50) 



l^nllc*[0,l] 



< cwu, 



n||C[0,l]) 



||C/n|]c[0,l] < C. 

Without loss of generality, there exists a limit s* G [0, 1] of the 
sequence {s n }. Assume that s* ^ s*, or, more particularly, s* > s* 
(the opposite case may be examined in a similar way). This yields 

that the value of the integral J 

o 

integration of implies that the left derivative (U, 



x) tends to zero as n — > +oo. But 
s n ) is equal to 
oo. 



/ 4> n (x)U n (x) dx. Due to f )50|) . this integral goes to zero as n — > 
o 

Using the information on the behaviour of U n summarized after equality 
( l4"6j) . we conclude that ||(C/n)a;||c[o,s n ] 0. Due to ( |50|) . without loss of 
generality there exists a constant U* such that 

(51) max U n (x) — > U*. 

xe[o,i] 



18 D.VOROTNIKOV 

But 

(52) max UJx) - min UJx) = U n (s n ) - UJO) 

xg[0,l] x6[0,l] 

< s n ||(C/„) a; || f7 [o iSn ] -»■ 0. 

Hence, 

(53) UJx) ->• £/* 

uniformly on [0,1]. 

Let us now test (Q with a smooth function h such that /i(s*) = 0, 
h(s*) = 1, h x (0) = 0, /^(l) = 0, obtaining 

i 

(54) J <f) n (x)U n (x)h(x) - U Tl (x)h xx (x) dx = h(s n ). 
o 

Passing to the limit, we find that the left-hand side goes to zero and 
the right-hand side goes to one, arriving at a contradiction. 

Appendix C. Proof of Theorem 14.11 
Consider the set of functions 

i 

]'{.r)/j(,r) d.r = I . (/'//},. > ,U. 

o 



B= {y(x)eL 1 (0,V 



y(x) > y(x + l/k) + My, for a.a. 0<x<l-l/k. 

Inverse induction shows that for any % — 1, . . . , k and y G B one has 

(55) N; > M. 

Let us define a mapping A on B. For each y G B, we let A(y) = Y, 
where Y is the solution of the problem 



k 

uY-qY xx = EW«„ XG(0,1) 
Y x (0) = Y x (l)~= 0. 



(56) 

To put it differently, 

k 

(57) Y = H 



i=i 
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Then, the set B is invariant for the map A. In fact, let y G B. Then 
(J55|) implies 

1 u 1 



J v(x)Y(x) dx = J2 ( v v)i = J 

n 1=1 n 



u(x)y(x) dx = 1. 



Further, 

i=\ i=l 

Finally, fix G [0,1 — 1/A;]. Then there is a number n such that 
x* G [x n , i ) . Set 



(58) JV* 



We claim that 

(59) ((^X - iV*) fe(x*) - + 1/A;)] > 0, z = 1, . . . , fc. 
Indeed, integration of ([22]) gives 

(60) ^(z) = / K 2 )?*^) <te, x<a u 



(61) Cft-cOz) = y v{z)qi(z)dz, x > a*. 

l 

Thus, the function g 4 is increasing on the segment [0, ai\ and decreasing 
on [ai, 1]. Assume first x* < a n . Then, if i < n, we have qi(x*) > 
qi(x* + l/k). Since y G 5, we also have > (^2/) n = -W*, and (1591) 

holds true. If i > n, we have qi{x#) < q^x* + l/k) and [yy) i < N*, and 
fl59|) again holds. If z = n, (159]) is trivial. Now, let x* > a n . In this 

case, if i < n+ 1, g;(x*) > ^(x* + and > {yy) n+l = iV*; if 

i > n + 1, qi(x*) < g«(x* + 1/A;) and [yy) i < AT*; and if i — n+ 1, ( 1591) 
is again trivial. 

Formulas (157]), (1551), (I26|) and (1551) yield 

fc 

K(x*) - Y(x* + l/k) = 22 ( u y)i[(li( x *) ~ Qi( x * + 

i=i 

> E - <li(x* + 1 A)] = N*[q(x*) - g(x* + l/k)} > M 7 , 

i=l 
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so the invariance of B is confirmed. 

Observe that A is a compact linear operator in Li(0, 1). Indeed, let 
B be the unit ball of the space Li(0, 1). Due to (jSTj) . its image A(B) 
is a bounded subset of the linear span of {<7i, • • • , <Zfc}, thus being a 
relatively compact subset of a finite-dimensional subspace of Li(0, 1). 

Let us show that (125]) may have at most one solution, so A can have 
at most one fixed point in B. If not, let Q be the difference of two 
distinct solutions. Then 

i 

(62) J u(z)Q{z)dz = 0. 

o 

Moreover, 

k ^ 

(63) Q = X>^> 

i=l 

whence 

fc 

(64) iyQ). = iyQ\iyS)p 3 = 1, • • • ,k. 

i=l 

From (1221) we deduce 



(65) / v(z) qi (z) dz = 1. 

o 

Therefore, the matrix [T 3 ^-] = (^ft),- is ergodic, i.e. it has positive 

entries, and the sum of the elements in every row is equal to one. By 
the Perron-Frobenius theorem, it has an eigenvector corresponding 

k 

to the simple eigenvalue 1, so that £j = an d & U ^ ne components 



=i 



is another 



£i are positive. On the other hand, by (|64|) . [Sj] = 

eigenvector of [Py] corresponding to the same eigenvalue. The sum of 
its components is zero due to (1621 . so it cannot be collinear with 

unless it is a zero vector. Since 1 is a simple eigenvalue, all {vQ) i are 
zeros, so Q = by virtue of (1631) . 

The set S is closed, convex and bounded in Za(0, 1). By Schauder's 
fixed point principle, A has a fixed point Q in 5, which is automatically 
a solution to (1251) . It remains to notice that Q is continuous as a linear 
combination of so (1271) holds for all < x < 1 — 
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Appendix D. Proof of Lemma 14.11 
Let us notice that 
(66) q = u + v, 

where u and v are the (unique) solutions to the following problems 



(67) 



(6* 



k 

vu - qu xx = J2$Si, xe(0, 1), 
u x (0) = u x (l) = 0, 

k 

vv - q,v xx = J2(Sai - 8 Si), xe(0,l), 
v x {0) = v x {l)~= 0, 



and 

(69) Si = S + Xi, i — l,...,k. 

Since S is the ///^-diffusive mean of and 1/k, there exists a solution 
U\ to the problem 

{vJJi - qU lxx = 8 Sl , xe(0,l/k), 
U lx (0) = U lx (l/k) = 0, 
U 1 (0) = U 1 (l/k). 

The solution u to flBTl) can be constructed in the following way: 

(71) u(x) — Ui(x — x^, Xi < x < Xi + 1/k, i — 1, . . . , k. 
Thus, u is 1/fc-periodic, i.e. 

(72) u{x) - u(x + l/k) =0, < x < x k . 
Set 

(73) d(x) = v(x) — v(x + 1/k), < x < x k - 
Then it suffices to show that 

(74) 7 = min d(x) > 0. 

Note that 

(75) vd - qd xx = 0, x G (0, x k ). 

By the maximum principle, if the minimum of d is non-positive, it is 
attained at or i^. To ascertain that this cannot happen, we are going 
to prove that 

(76) 4(0) < 0, d x (x k ) > 0. 
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Set 



X 



(77) V(x) = / v(z)v(z) dz, 







and let # be the solution of the Cauchy problem 



k 



0* = 52(^-5%), 16(0,1), 



i=l 



0(0) = 0. 



Note that 9 is non-negative and l/fc-periodic. 
Integration of (|68l) gives 



(79) ^ - Wb = 0. 
Therefore 

(80) 4 = "(*>-^ + V*>. 

From (179")) we deduce 

(81) V(0) = V(1) = 0, 
and 

(82) v-q(-) =9>0. 



Using the distributional maximum principle [TBI Theorem B], we con- 
clude that 

(83) V{x) > 0, < x < 1, 



so 



(84) 4(0) = -W < o, 4W = > 0. 



Appendix E. Proof of Lemma 15.11 
It suffices to prove that the set 

O = ^ + ^|^6C 2 [0,l],^(0) = ^(l)=0} 

is dense in C[0, 1]. 

Let h G C 2 [0, 1] be an arbitrary function which is locally constant 
near the zeros of b. These functions constitute a dense subset 0\ of 
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C[0, 1]. Let 



y 

(85) </?(x) = h(x) + / exp I / dt I h y (y)dy, 



Kbit) 



Clearly, is equal to a constant c i (resp. c^~) in a left (resp. right) 
neighbourhood of the point X{. But 

(86) h = ip + Kbip x , 

so c~ = cf = h(xi). Thus, <p is C7 2 -smooth and y9 x (0) = y? x (l) = 0. By 
virtue of (1561) . Oi is contained in O. 



Appendix F. Proof of Lemma 15.21 
The solution r to (!30|) can be written explicitly: 



(87) r(x) 



X / X \ 

■— — — / exp / — — dt R(s) ds, Xi < x < a. 



r ( x ) = n \ I ex P I / -rzrv ^ I -^l 5 ) ds, di< x < x i+ i. 
k>K x ) J \J ] 



Note that it is unique in Li(0, 1). Indeed, if r\ is another solution, then 
w = r — ri satisfies the conditions of Lemma 15.11 

We need to show that r — >■ ^ -Rj5 ai weakly-*, uniformly with respect 

i=l 

to i?. It suffices to prove that, for each i, r — > Ri5 ai weakly-* on the 
interval (xi,x i+ i), uniformly in R. We restrict ourselves to the case 
i = 1, and the others are analogous. 
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We calculate, integrating by parts, 



l/k 

(89) r 1 = J r(x) 



dx 



a x 1 x \ 





l/k l/k 
+ I ; / \ / exp 



nb(x) 



Kbit) 



dt R(s) ds dx 



a x 



x / x 



-\ 



/ exp / JtT) dt ) R{s)ds 



+ 



R(x) dx 



o 



/ exp / 



l/k 



Kb{t) 



dt R{s) ds 



i/k 



l/k 



l/k 



\s 



Let us show that for every x* e (0, a) 



(90) 



lim / r(x) dx = 0, 

re— >+oo y 



uniformly in i?. Indeed, let s K < rr* be such that 



(91) 



Sw) dt=v ~ K - 
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(92) / r(x)dx 



X* XIX \ 

/^)/ exp / ^T) dt ) R{s)dsdx 



1 



cxp 



1 



< 



1 — exp 



Kbit) 



Kb{t) 
1 — exp 

dt 



dt R(s) ds 



dt 



+ / R(x)dx 



nb{t) 



R{s) ds+ 



R(s) ds 



1 — exp 



nb{t) 



dt 



< J R{s) ds + [1 - exp(/sT 1/2 )] j R{s) ds -> 



as k — > +00. 

Due to (ED]), for any / G C[0, 1/Jfe], 



(93) 



.t'» 

lim / r(x)fo(x) dx = 0, 



uniformly in R. 

Similarly, for all x* G (a, 1/fc) and / G C[0, 1/fc], 



(94) 



i/fc 

lim / r(x)fo(x)v = 0, 



uniformly in R. 

Fix e > and / G C[0, l/k]. Let x* and x* be so close to a that 
\f(x) — /(a) I < e/2 provided x* < a; < x*. Then 



(95) 



r(x)[f(x) - f(a)} dx 



< e/2. 
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Due to (E3D and (El with / = / - f{a) 



(96) 



1/k 



r(x)[f(x) - f(a)]dx- / r(x)[f(x) - f(a)] dx 



< e/2 



for sufficiently large k. Thus, 



i/fc 



(97) (r- J R 1 «J 01 ,/) = (r-fi«J a ,/) = y r(x)[/(x) - /(a)] < e. 

o 

Appendix G. Proof of Lemma 15.31 
Let r aK be the solution of the system 



(9£ 



- «(&»V,*)a> = i^ 1#t , X G (0, 1) 

r CT , K > 0. 



Then, by Lemma [5.21 

(99) lim d(r ffl(( T(S3 j i ai ) = ) 



1=1 



uniformly in cr. Thus, it suffices to prove that for every k there is e K > 
such that 

lim d(w a>K ,r aK ) = 0. 

K— >+oo, cr<e K 

This would follow from the claim that for every k there is e K > so 
that for a < e K we have d(w,j, K , r ffjK ) < 1/k. If it is not true, then for 
some k there exists a sequence a n — > such that 

Since w fii)K and r CTn A are solutions of the problems (!32j) and (1981) . we 
have 

(100) - a n (w an , K , (p xx /r]) + (w UntK , y? + nb(p x ) = (R an , K , <f), 

(101) (r antK , if + Khp x ) = (R an ,K, <p), 

for any cp G C 2 [0,1], <p x (0) = <p x {l) = 0. Since the sequences w anyK 
and r a K lie in the space of probability measures, which is weakly-* 
compact, without loss of generality there exist their weak-* limits w K 
and r K . Clearly, 



(102) 



d(w K ,r K ) > 1/k. 
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On the other hand, taking the difference of f llOOp and (llOip . and passing 
to the limit, we find (w K — r K , if + Kb(p x ) = 0, so w K = r K by Lemma 
5.11 and we arrive at a contradiction. 



Appendix H. Proof of Lemma 15.41 

The pair (w a , K , R^k) = (t)P, vP) satisfies ( 1521) . so, by Lemma 15751 for 
every k there exists e K > such that 

k 

(103) Jim d( W ,X>^)A;) = 0- 



K— >+oo, a<t K 

i=l 



Multiplying the second equation in (1T51) by P and integrating, we 
find 

l l 



(104) / v(x)P\x) - qP xx (x)P(x) dx = \ rj(x)p(x)P(x) dx, 

whence 



(105) inf u(x) / P 2 (x)dx + q / P£(x)dx< sup F(ac). 

0<x<l J J 0<x<l 

o o 

Hence, 

(106) 11^11^(0,1) < C||P||qp,i] < C\\P\\ w i m < C. 

Assume that (1541) is not true, i.e. there exist 5 > and sequences 
K n — > oo and <r n < e Kn such that for the corresponding solutions 

(p„,P„) = (P(r nj « B ,-Per n ,K„) to ^TSJ) we have ||P„-Q||c[o,i] > Since tne 
embedding W^O, 1) C C[0, 1] is compact, without loss of generality we 
may assume that P n converges to some limit P in C[0, 1]. Obviously, 

(107) ||Po-Q||c[o,i]><5. 

Passing to the limit in the second, forth and the last equations in (1T6|) 
- the combination of the first two is understood in the weak sense (JBJ) 
— and remembering ( I103p . we find 

k 

i=l 

(Po)x(0) = (Po)x(l) = 0, 
1 

i/(x)Pq(x) dx = 1. 
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By Theorem 14.11 P coincides with Q, which contradicts (II 07ft . 
From (I103p we deduce 

(108) lim d (p, y ^1^-) = 0. 

Due to and l//c-periodicity of r), (11081) implies that 



k 



(»Q)A, 



(109) ,(«) ' 

weakly-* as k — > +oo, cr < e K . Taking test functions which are equal 
to 1 in one of the wells and are zero at the minima of the potential 
located outside of that well, we derive (1351) from (1 1 9 H . 
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